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ABSTRACT 

We discuss a new one-dimensional non-LTE time-dependent radiative-transfer technique for 
the simulation of supernova (SN) spectra and light curves. Starting from a hydrodynamical 
input characterizing the homologously-expanding ejecta at a chosen post-explosion time, we 
model the evolution of the entire ejecta, including gas and radiation. The boundary constraints 
for this time-, frequency-, space-, and angle-dependent problem are the adopted initial ejecta, 
0^ a zero-flux inner boundary and a free-streaming outer boundary. This relaxes the often un- 

. suitable assumption of a diffusive inner boundary, but will also allow for a smooth transition 

from photospheric to nebular conditions. Non-LTE, which holds in all regions at and above the 
photosphere, is accounted for. The effects of line blanketing on the radiation field are explic- 
itly included, using complex model atoms and solving for all ion level populations appearing 
in the statistical-equilibrium equations. Here, we present results for SN1987A, evolving the 
model "lml8a7Ad" of Woosley from 0.27 to 20.8 d. The fastest evolution occurs prior to day 
1, with a spectral energy distribution peaking in the range ^300-2000A, subject to line blan- 
keting from highly ionized metal and CNO species. After day 1, our synthetic multi-band 
light curve and spectra reproduce the observations to within 10-20% in flux in the optical, 
with a greater mismatch for the faint UV flux. We do not encounter any of the former dis- 
crepancies associated with the He I and H I lines in the optical, which can be fitted well with 
a standard Blue-supergiant-star surface composition and no contribution from radioactive de- 
cay. The effects of time dependence on the ionization structure, discussed in Dessart & Hillier, 
are recovered, and thus nicely integrated in this new scheme. Despite the ID nature of our ap- 
proach, its high physical consistency and accuracy will allow reliable inferences to be made 
CO ' on explosion properties and pre-SN star evolution. 
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, 1 INTRODUCTION SN light, using photometric, spectroscopic, or spectropolarimetric 
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data. Hence, developing accurate radiative-transfer tools capturing 
the key physics controlling the interaction of light and matter is of 
prime importance. 



Core-collapse supernovae (SNe) are extraordinary events situated 
at the crossroads of many fields of astrophysics. They mark the 
birth of compact objects, either neutron stars or stellar-mass black 

holes which, owing to their high compactness, are often the site Of all SNe, only core-collapse, and in particular Type II SNe, 

of tremendous magnetic fields, as in magnetars, or the site of have well-identified progenitors. Numerous Type Il-Plateau (II-P) 

tremendous rotation rates, as in millisecond-period pulsars. Their SNe have now been associ ated with the explosion of low-mass 

ejecta make a significant contribution to the chemical enrichment Red-Supergiant (RSG) stars ( Smartt 2009), while the progenitor of 

of galaxies, while affecting their dynamics and energetics. The con- the Type II-peculiar SN1987A was a Blu e-Supergiant (BSG ) star, 

nection to 7-ray bursts, for a subset of these SNe, may turn them named Sk -69 202 (for a review, see, e.g., lArnett et alj|l989h . The 

into excellent probes of the early Universe. Our focus, however, excellent quality of the observational data for SN1987A makes it 

is to characterize the SN ejecta itself and extract information that ideal for detailed modeling. Numerous radiative-transfer studies of 

can help us understand the properties of the explosion and of the the e arly-time spectra of SN1987A have been done, with some suc- 

progenitor star. Such inferences are based on the analysis of the cess (Eastma n & KirshneJl98 9; Hoflich 1987l;lLucvll987j: | Hoflichl 

[l988! - |Luc-vlll987l : ISchmutz e7alj|l99(ll : lMazzali et alj|l992n . How- 
ever, these studies are now twenty years old. They generally as- 

* E-mail: Luc.Dessart@oamp.fr sumed Local Thermodynamic Equilibrium (LTE; at best treating 
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only a fraction of the species/levels in non-LTE) and/or steady- 
state for_the radia tive transfer, and em ployed small model atoms. 
iLucvl il987t) and ISchmutz et alJ dl99Clh have emphasized the ef- 
fects of line blanketing in the UV. One explicit question raised by 
these studies was the problematic observation of He I lines in opti- 
cal spect ra, which required unaccepta ble helium enrichments. More 
recently, iMitchell et all ( 1200 ll l2002h argued for significant mixing 
of 56 Ni beyond 5000 km s _1 in order to reproduce Balmer line pro- 
files during the first weeks after explosion. 

The Balmer line strength problem has now been associated 
with the erroneous neglect of important time-dependent terms 
that appear in the energy and statistical-equil ibrium equations 
dUtrobin & Chugail2005l ; IPessart & Hilliej2 008), while reproduc- 
ing the characteristics of He I lines seems to require a non-LTE 
treatment. Hence, some of these discrepancies may simply reflect 
the shortcomings of the radiative-transfer tools employed rather 
than a genuine peculiarity of the progenitor. Taking a new look 
at this dataset therefore seems warranted. Furthermore, studying 
SN1987A is a good exercise to gauge the level of accuracy of 
radiative-transfer codes and assumptions, facilitated by the consid- 
erable advances in computer technology over the last twenty years. 
Because atomic data represent a fundamental and essential ingre- 
dient of radiative-transfer computations, the considerable improve- 
ments in that domain over that period make such calculations more 
accurate than calculations done when SN1987A went off. 

In this paper, we present a new approach for non-LTE time- 
dependent radiative transfer modeling of SN ejecta using the cod e 
CMFGEN dHillier & Miilerll 19981 : IPessart & Hillie"ill2005bl I2008T) . 
We discuss the conceptual aspects of the method and emphasize 
the key equations that are solved, delaying a more comprehensive 
presentation of the technical details to a forthcoming paper (Hillier 
& Dessart 2010, in preparation). We illustrate this new capabil- 
ity with results obtained for SN1987A, starting our time evolution 
from an hydrodynamical input of an exploded BSG star (Woosley, 
priv. comm.; model "lml8a7Ad"). We do not present an in-depth 
study of SN1987A, but merely use this well-observed SN to check 
and confron t our model results. A p reliminary set of results were 
presented in IPessart & Hillier! (2009). In the following section, we 
discuss the various approaches we have used in the recent years for 
SN spectroscopic modeling, emphasizing their merits and limita- 
tions. We then describe in £[3]the hydrodynamical model we employ 
as a basis for this calculation, summarizing its various properties, 
before presenting our setup for the radiative-transfer calculation. In 
Sj4] we describe the ejecta evolution, such as temperature and ion- 
ization structure, and present checks on our numerical technique. 
We then discuss in Sj5]the radiative properties of our non-LTE time- 
dependent models, covering in turn synthetic spectra and the bolo- 
metric light curve, and detailing in particular the sources of line 
blanketing at various epochs. In ij6] we compare our theoretical 
predictions to UV and optical observations of SN1987A. Finally 
we present our conclusions in S|7] and lay out the various projects 
ahead. 



2 THE RADIATIVE TRANSFER PROBLEM AND THE 
APPROACHES TO ITS SOLUTION 

Radiative-transfer modeling of SN ejecta can be performed in a va- 
riety of ways and with different accuracies and consistencies: the 
interaction between the gas and the radiation may be done in LTE 
or in non-LTE; the problem may be simplified to its steady-state 
form (c — ¥ oo and d/dt =0); the relative importance of absorption 



versus scattering opacity may be imposed or calculated. The latter 
is controlled by transition rates, level populations, and the ioniza- 
tion state, which are all connected to the radiation field and thus the 
whole problem is tightly coupled. 

Our fundamental aim is to model the time evolution of the gas 
and radiation field from near the onset of the explosion (i.e., shock 
breakout in core-collapse SNe) when the ejecta are hot, dense, and 
possess a large optical depth, until the nebular phase when expan- 
sion and cooling have made the ejecta cool, tenuous, and optically 
thin at all wavelengths. To do so, we need to track the evolution of 
the following: 

1) The specific intensity 7 as a function of radius, time, fre- 
quency, and angle, i.e., I = I(r, t, v, fj,), as well as its moments J, 
H, and K. 

2) The properties of the gas, such as electron density n c and 
ion density ni, temperature, energy. This requires accounting for 
changes in excitation and ionization, as well as for contributions 
from radioactive decay of unstable isotopes. 

3) The coupling between the radiation and the gas, i.e., opac- 
ity and emissivity. This represents a considerable task in SN ejecta 
because of departures from LTE conditions and the importance of 
scattering. 

4) The changes in composition of unstable isotopes. 
Finally, all this work has to be done using a reliable and compre- 
hensive atomic dataset for a wide range of species and ions. 

Schematically, assuming homologous expansion and retaining 
only the terms of order v/c, the combined set of equations we wish 
to solve are: 

• The statistical equilibrium equations (one for each level i of 
each ion treated): 

Dm/p 1 D(r s m) ^ I n u ^ , n 

j 

• The energy equation: 

De P Dp f ge dccay 

~Dt ~ ~~p~Dt ~ J MXvJv ~ ^ + ; (2) 

• The time dependent zeroth and first moment of the radiative 
transfer equation: 

1 D{r 3 J v ) 1 d(r 2 H„) vV dJ y 



cr 3 Dt 
and 
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= Vv — XvJv (3) 
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(4) 

-XvHv . (5) 
rc av 

In the above equations D/Dt is the Lagrangian derivative, rii the 
number density of state i, Rij is a term that represents the rates 
connecting the level i to the level j, e is the internal energy per unit 
mass, v the frequency, \v the opacity, r\ v the emissivity, J v the 
mean intensity, and edccay is the specific energy associated with 
the decay of unstable nuclei. The internal energy can be written in 
the form 



e = ex + ei , 
where 

3kT(n + n e 



en 



2pmn 



(6) 



(7) 
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ei 



limn 



(8) 



n is the total particle density (excluding electrons), m is the atomic 
mass unit, fi is the mean atomic weight, and Ei is the total energy 
(excitation and ionization) of state i. 

In a steady-state approach, we drop all terms involv- 
ing partial-derivatives with respect to time, in the statistical- 
equilibrium/energy-equations and/or in the moments of the trans- 
fer equation. By neglecting these in both, the approach allows for 
departures from LTE and assumes steady-state ( £|2. 1 1 . Retaining 
them in the first but neglecting them in the second, one grasps the 
time-dependent effects important for the ionization structure, while 
the radiation field is still assumed steady-state (*j2j2j- By retaining 
these terms in all equations above, one accounts for the full time de- 
pendence of the problem and thereby reaches a much higher level of 
physical consistency ( i]2,3t . These various approaches correspond 
chronologically to what we have followed so far. We now present 
them in more details. 



2.1 Non-LTE and steady-state approach 

Because Type II-Plateau (II-P) SNe, associated with the explo- 
sion of RSGs, are characterized by a massive quasi-homogeneous 
hydrogen-rich homologously-expanding ejectum during their pho- 
tospheric phase, their photospheric properties are weakly sensitive 
to the details of the explosion mechanism or morphology at early 
times. In particular, the hydrogen-rich region of Type II-P SN ejecta 
seems to remain quasi-spherical even when there is strong evidence 
for an aspherical explosion, as revealed by the growing linear po- 
larization of continu um photons at the end of the plateau phase 
dLeonard et alj|2006l) . This property gives impetus to model atmo- 
sphere calculations in which multi-dimensional effects are ignored, 
but which include a much more accurate treatment of non-LTE 
and line-blanketing effects (both of which are very computationally 
costly). Being free of these various uncertainties, Type II-P SNe are 
ideal for benchmarking purposes, and we therefore focused our ini- 
tial investigations on these objects. 

Our initial perspective was to consider the photospheric lay- 
ers of the SN ejecta to be just like those of a stellar atmosphere, 
characterized here by a homologous and fast expansion, a steep 
density distribution, and spherical symmetry. In this approach, 
we neglected the Lagrangian derivative terms in the radiative- 
transfer, statistical-equilibrium, and energ y equations and solved 
the radiative-transfer pr oblem in ID (see Hillier & Milleri [l998l : 
iDessart & Hillier 2005b)Q Observ ed spectra are computed in the 
observer's frame using the code of iBusche & HillieJ 12005), with 
relativistic terms taken into account but ignoring time dependence 
(see Hillier & Dessart 2010 for details). 

The radiative-transfer grid is set to cover the vicinity of the 



1 If we neglect only dJ/dt we get extra terms of order v/c in the transfer 
equation which is generically referred to as the relativistic form of the trans- 
fer equation. These extra terms are unimportant in O and Wolf-Rayet stars, 
but are of crucial importanc e for the temperature stru cture in SN, particu- 
larly below the photosphere (Mihalas & Mihalas 1984). Various approxima- 
tions to the transfer equation have been suggested to avoid the explicit-time 
dependence in the trans fer equation, but these are ultimately unphysical. 
Pinto & Eastman 12000) concluded "that a fully time-dependent solution to 
the transport problem is needed in order to compute SNe la light curves and 
spectra accurate enough to distinguish subtle differences of various explo- 
sion models", and we believe the same must be true for Type II SNe. 



photospheric region, extending inwards to a Rosseland mean opti- 
cal depth of ~100 and extending outwards to distances where ra- 
diation is in a free- streaming regime. In between the boundaries, 
we adopt an analytical distribution for the SN ejecta density. We 
impose a uniform composition (i.e., the modeled ejecta layers are 
chemically homogeneous), so that changes in ion densities reflect 
exclusively changes in ionization. The temperature structure of the 
ejecta is controlled through the imposed diffusive inner boundary 
flux and constrained through fits to the spectral energy distribution 
(SED) as well as line profile strength and shape. An independently 
inferred distance to the Galaxy host and/or the estimated time of 
explosion allow one to approximate the radius at a given epoch, 
thereby fully constraining the problem. 

Equipped with the assumptions outlined above, we analyzed 
the photospheric phase of SNe 1999em, 2005cs, and 2006bp 
( IDessart & HillieJ (20061 ; iBrown et all l200l IDessart et alj l2008h . 
We discussed the importance of line blanketing on the SED, the im- 
portance of non-LTE effects for line formation, and the importance 
of line overlap for proper identifications of spectral features. We 
found that at early times, prior to the ejecta recombining to a neu- 
tral state, our approximations produced very satisfactory fits to ob- 
servations, with ejecta properties compatible with the RSG nature 
of the progenitor star. Based on these quantitative analyzes, we also 
determined distances to these SNe, obtaining values in agreement 
with those obtained with independent methods (e.g., Cepheids). 

The advantage of this method is that it is self-contained, so that 
one can compute a model without considering the previous evolu- 
tion or the ejecta conditions below the inner grid boundary. But 
this approach becomes inaccurate when the time-dependent terms 
appearing in the equations above become non-negligible, or when 
the inner boundary flux cannot be accurately described by a steady- 
state diffusion solution. 



2.2 Non-LTE radiative transfer with steady-state 

radiation-field but time-dependent terms in the rate and 
energy equations 

The steady-state assumption is a simplification of the radiative 
transfer problem, one that holds safely in the context of hot star 
atmospheres for which CMFGEN was originally designed. How- 
ever, in SNe, low density and fast expansion conspire to make 
the expansion and the recombination timescales eventually com- 
parable. Exacerbated further by optical-depth effects, time de- 
pendence can lead to an ejecta ionization freeze-out, modifying 
the electron and ion densities, line and continuum optical depths, 
and conseq uently the emergent flux globally. This was demon- 
strated by lutrobin & Chugail fe005l) in the case of SN 1987A, 
with primary focus on Ha and BaIl6142A. By retaining the time- 
dependent terms in the statistical equilibrium and energy equa- 
tions in CMFGEN (while keeping the same approach otherwise), 
we confirmed and generalized the strong time-dependent effects on 
Type II SN ejecta ionization, and extended the i mplications to the 
entire spectrum, including continuum and lines ( Dessart & Hillier 
2008). In particular, we demonstrated the importance of these time- 
dependent terms for reproducing the strength and shape of Balmer 
line profiles, although we emphasized that all lines and all ions are 
to some extent affected. Even at early times when hydrogen is fully 
ionized, time-dependent effects on the helium ionization are no- 
ticeable, leading to stronger and broader He I lines compared to 
equivalent steady-state models. 

While such impor tant time-dependent terms were included in 
Dessart & Hillier! d2008l) . the approach focused only on a portion 
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of the SN ejecta, characterized by a prescribed density distribution 
and chemical homogeneity, with an imposed flux at the inner, diffu- 
sive, boundary. While suitable for general studies of the radiative- 
transfer problem, and in particular for gauging the importance of 
time dependence on the ejecta ionization structure, this setup is not 
very amenable for comparing models with observations, or for as- 
sessing the adequacy of a given explosion model. 

2.3 Non-LTE approach with full time dependence 

The recent improvements to the code retain all the assets of the 
previous version, namely non-LTE and time dependence in the 
statistical-equilibrium and energy equations. However, in addition 
to the derivatives with respect to space, angle, and frequency, we 
now treat time derivatives of the moments of the mean intensity 
that control the transport of radiation through the ejecta. For sim- 
plicity, we still ignore time derivatives of the specific intensity. Im- 
portantly, rather than focusing merely on a portion of the ejecta, we 
start from hydrodynamical inputs describing a given explosion and 
we map the entire ejecta onto the radiative-transfer grid. We then 
follow explicitly the evolution of both the radiation and the gas at 
all depths. Hence, with such a "full-ejecta" approach, we now re- 
tain all time-, depth-, angle-, and frequency-dependent terms in the 
radiative-transfer equations, the energy equation, and the statistical- 
equilibrium equations. Importantly, this approach will allow us to 
model the full non-LTE spectral evolution as a function of time, 
which in turn will allow us to accurately model the flux and color 
evolution of the hydrodynamical model of the SN. 

The treatment of boundaries is more straightforward in this 
approach. The outer boundary is positioned at large distances (cor- 
responding to a maximum velocity as high as ~60000km s _1 ), 
where the radiation is freely streaming at all wavelengths apart 
from, perhaps, the He I-He II continua - this matters little since few 
(if any) photons are injected at such wavelengths at any time. The 
inner boundary flux is constrained by imposing a zero-flux condi- 
tion. This is a reasonable approximation when the entire envelope 
is optically thick — it ignores effects due to Doppler shifts at the 
inner boundary and the finite light travel time. At later times we 
use a a nebular approximation (the incoming and outgoing specific 
intensities are equal, modulo Doppler shifts which we correct for). 
Diffusive energy emerging from large optical depths as time pro- 
ceeds is naturally taken into account. The parameters describing 
the model are entirely determined at the first time step, and we just 
let the ejecta evolve under the various conservation laws for level- 
population rates, charge and gas/radiation energy. This is primar- 
ily an initial- value problem, with straightforward boundary condi- 
tions at spatial boundaries. Full details of the solution technique, 
and approximations, will be provided by Hillier & Dessart (2010, 
in preparation). 

Using a hydrodynamical model as input for our radiative- 
transfer calculations represents a big improvement over our for- 
mer approaches, which required the setting of all ejecta conditions 
(density, radius, velocity and boundary conditions). The past free- 
dom allowed the tuning of model parameters but lacked physical 
consistency. Although it will now be harder to fit a given observa- 
tion, this new approach allows a more direct confrontation of syn- 
thetic observables with ejecta and progenitor properties, providing 
constraints on the progenitor star and its pre-SN evolution, and the 
characteristics of the explosion. 

To track composition changes with depth, we allow for a 
depth-dependent composition. In addition, we follow the change 
in composition and the associated energy release following the de- 
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Figure 1. Elemental mass fractions as a function of Lagrangian mass for the 
hydrodynamical input model "lml8a7Ad" used as the starting ejecta con- 
ditions for our radiative-transfer study. We show hydrogen (black), helium 
(blue), oxygen (orange), carbon (green), and 56 Ni (red). We also overplot 
the velocity as a dotted black line. Note that there is ~ 1 . 1 Mq of ejecta 
material traveling faster than 4500 km s — which is the photospheric ve- 
locity at the end of our time-sequence at 20.8 d. The whole spectroscopic 
and photometric evolution presented in this paper thus corresponds to a con- 
stant photospheric composition, i.e., that of the progenitor surface. 



cay of unstable nuclei. Since most explosion models give J Ni as 
the most significant radioactive species, we follow only one decay 
chain. We assume either pure local energy deposition or solve for 
the non-local deposition using a Monte-Carlo approach for 7-ray 
transport (this code will be presented in a separate paper). In the 
present calculations on the early-time evolution of SN 1987A, the 
typical 7-ray photon mean-free path is much smaller than the scale 
of the ejecta so that the energy deposition is purely local. 

In the following section, we illustrate the results from our non- 
LTE time-dependent radiative-transfer modeling for SN1987A, 
based on the model "lml8a7Ad" (Woosley, priv. comm.). With our 
approach, we can now compute both spectra and light curves in 
full non-LTE and allow for time dependence, with explicit account 
of line blanketing. Such an approach has never been followed in 
the past and constitutes a key asset of our work. In the present pa- 
per, we wish to highlight the key properties of the SN1987A ejecta 
and show the match of our synthetic spectra to those observed. 
Aiming at validating our code against the excellent observations 
of SN1987A, we do not address specifically the issues concern- 
ing the explosion, the ejecta geometry, or the progenitor identity 
of SN1987A. Such an in-depth study will be performed later, al- 
lowing for variations in the properties of the progenitor and of the 
explosion, and gauging the impact on observables. 



3 INITIAL HYDRODYNAMICAL MODEL AND SETUP 

The hydrodynamical model we use for the calculations presented 
in this paper was produced in two steps by Stan Wo osley (priv. 
comm.), using the code KEPLER dWeaver et alj 1 19781) . First, an 
18 Mq star on the main sequence was evolved with allowance for 
rotation, magnetic torques, and mass loss. The angular momentum 
halfway through hydrogen-core burning on the main sequence is 
2.9 x 10 52 erg-s. At the stellar surface, the rotation speed at that 
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Table 1. Ejecta composition, given in solar masses, of the "lml8a7Ad" 
model used as input for our radiative-transfer computations. Note that all 
metal mass-fractions have a floor value corresponding to the LMC metal- 
licity. The total ejecta mass is approximately 15.6 Mq 



Element 


Mass [M Q ] 


Element 


Mass [M©] 


H 


3.94 


Si 


0.08 


He 


9.04 


S 


0.03 


C 


0.61 


Ar 


0.005 


N 


0.027 


Ca 


0.004 


O 


1.2 


Ti 


6.5xl0- 5 


Ne 


0.34 


Fe 


0.01 


Mg 


0.06 


56 Ni 


0.084 



time is 240 km s" 1 on the equator, which corresponds to ~33% of 
the local Keplerian (critical) velocity. An LMC metallicity Z was 
adopted, with Z = OAZq. No convective overshooting was used, 
but a modest semiconvection was allowed for. 

At the onset of collapse, the star is a blue-supergiant 
(BSG) with a radius of 3.28 xlO 12 cm (~47R©), a luminosity of 
8x 10 38 erg (~2.1x 10 8 L Q ), an effective temperature of 18000K, 
and has an age of 14.7 Myr. 

The pre-SN envelope above a mass cut of 1.47 M© (the 
proto-neutron star which forms a pulsar with a rotational period 
of ~12ms) is then exploded with a piston. Its trajectory is char- 
acterized so that it yields an ejecta mass of 15.45 M© with an 
as ymptotic kinetic energy o f 1.2xl0 51 erg (for details, see §2.3 
of IWooslev & Weaver! fl995h . The ejecta composition we inherit, 
when we start our radiative-transfer computations, is given in Ta- 
ble Q] while the distribution of the species with mass and velocity 
is shown in Fig. Q] The model provided did not include a com- 
prehensive description of isotopes, and lacked some elements with 
low abundances (e.g., sodium). For these, we adopted the corre- 
sponding LMC-metallicity value (given by scaling the solar metal- 
licity value by 0.4) and set it to be constant with depth (we ex- 
pect an abundance variation at the photosphere as it recedes in the 
ejecta towards deeper mass shells, but not over the ~21 days cov- 
ered by our time sequence so this expedient should be adequate). In 
practice, we use a depth-independent mass fraction of 1.35x 10 
for sodium. Similarly, although 44 Ti is present and quite abundant 
at depth in the ejecta, it is absent at the surface and so we en- 
force a floor value of 1.32x 10 -6 for titanium to reflect the LMC 
metallicity value for that species. Barium and Scandium, which are 
both strongly under-abundant, are ignored in this first set of sim- 
ulations, and hence their associated line features are absent from 
our synthetic spectra. Note that some chemical mixing has been 
done to mimic the effects of multi-dimensional fluid instabilities 
(e.g., Rayleigh-Taylor). The moderate mixing softens the composi- 
tion gradients but preserves the strong ejecta stratification, hydro- 
gen being absent in the inner 3 solar masses of the ejecta (below 
1500 km s^ 1 ). The 56 Ni mass fraction is sizable only below 1000- 
2000 km s _1 but owing to mixing, shows a non-zero mass fraction 
even at the outer edge of the ejecta. This is unphysical but tests we 
performed show that such a low 56 Ni abundance does not impact 
the gas, nor the radiation, properties in regions exterior to a velocity 
of 4500 km s _1 (the ejecta photospheric velocity at the end of our 
time sequence) during the 21 -day timespan we cover. 

Before starting our radiative-transfer calculations, we need to 
make a number of adjustments to this input hydrodynamical model. 
The innermost part of the ejecta is infalling at moderate speeds 



and shows a pronounced broad peak in both temperature and den- 
sity. We thus start our grid just outside of that region, at a radius 
of 1.3xl0 12 cm corresponding to a velocity cut of 500km s _1 . 
This excision, which amounts to ~0.02M©, is very small. As 
the outer regions of the ejecta are described by only a handful of 
points beyond ~ 10000 km s~ , and reach a maximum velocity of 
a third of the speed of light, we trim this poorly-resolved region 
at ~8000km s , and stitch an outer region in homologous ex- 
pansion, constant temperature and a power-law density distribu- 
tion with an exponent —20, extending to a maximum radius of 
1.3x 10 14 cm and a maximum velocity of 55000km s _1 . Although 
at the time of 0.27 d in this initial model homologous expansion 
is pretty well established (apart from the inner regions which we 
have trimmed, r/v varies by about 15% between the inner and 
the outer ejecta), we enforce strict homologous expansion (exact 
numerically, i.e. r /inconstant), using the values of velocity and 
radius halfway through the ejecta. This is done because our time- 
dependent radiative-transfer solver requires dv/dr — v/r, but this 
tinkering does not alter the corresponding time of explosion by 
more than a few percent in the surface layers we model here. That 
homology does not hold exactly in the ejecta at 0.27 days results 
primarily from the finite size of the progenitor star, whose sur- 
face radius represents 10% of the SN photosphere radius at 0.27 d. 
These edits to the initial model are merely for starting our time- 
sequence. After that, we just let the ejecta evolve under the con- 
straints set by radiation transport and energy/charge conservation. 

To limit the size of the numerical problem to solve, we ex- 
clude some species/ions that do not contribute sizably to the prob- 
lem, e.g., species that are not identified conspicuously in spectra 
at any time (e.g., neon or sulfur), or ionization stages that are not 
needed for the accuracy of the problem (e.g., high ionization stages 
of sodium, magnesium, or calcium). This expedient is justified at 
this exploratory stage, in particular given the brute-force approach 
followed. In the future, we will employ a more complete atomic 
model. In Table [2] we provide details of the model atom used in 
the present computations (the same one is used throughout the se- 
quence), including for each ion species the number of super and full 
levels, the number of transitions, and the uppermost level included. 
Overall, we identify only one inadequate setting in the present cal- 
culation, which has to do with the neglect of Ni II, Ni III, and Ni IV 
in the model atom. 

The basic a tomic processes, atomic data, and its utilization, 
are discussed by Iffillier & Millerl dl998h . In the code we include 
the following processes: bound-free, bound-bound, free-free, di- 
electronic recombination, electron scattering, Rayleigh scattering 
by hydrogen, inner-shell ionization by X-rays, collisional ioniza- 
tion/recombination, collisional excitation/de-excitation, and charge 
exchange reactions with H and He. To simplify the atomic model 
atoms we utilize super-levels, although all transitions are still 
treated at their correct wavelength. Photoionization cross sections 
are smoothed (typically with a Gaussian of width 3000 km s _1 ) 
with recombination treated as the inverse of photoionization. At 
present we do not include ionization by non-thermal electrons 
which are particularly important at later times. 

For most ions there exist states that lie above the ground state 
of the ion, and these can be treated in a variety of ways. Generally, 
for states permitted to autoionize in LS coupling, they are treated 
as reson ances in the p hotoionization cross-section. The Opacity 
Project dSeatonl Il987h data explicitly contain these resonances 
while older calculations may not. The treatment of states not 
permitted to autoionize in LS coupling varies - in some cases (e.g., 
C III) these states are expected to autoionize and are added to the 
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Rphot Vphot T phot n e phot 




Table 2. Summary of the model atom used in our radiative-transfer calcu- 
lations. The source of the atomic datasets is given in text. refers to the 
number of full levels, N s to the number of super levels, and Nt rans to the 
corresponding number of bound-bound transitions. The last column refers 
to the upper level of a given ion included in our treatment. At the bottom 
of the table, we give the total number of full levels treated, and the corre- 
sponding number of transitions explicitly included. 

Species N 3 Nt rans Upper Level 



1 10 
Time since Breakout [d] 



Figure 2. Evolution of the radius (black; Hm = R/10 15 cm), veloc- 
ity V 4 (blue; V 4 = V/10 4 km s" 1 ), temperature T 4 (red; T 4 = T/10 4 K), 
and electron density n e ,11 (green; n e ,n = rie/10 11 cm -3 ) at the photo- 
sphere as a function of the time since shock breakout (we adopt for the pho- 
tosphere location the point where the inward-integrated electron-scattering 
optical depth is 2/3). We use a logarithmic scale for the abscissa to bet- 
ter render the fast early evolution of the photospheric conditions. The dots 
correspond to the actual times at which the computations were performed. 



photoionization data. In other cases the levels are ignored although 
for many CNO atoms we now have accurate autoionizing rates 
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&A collaboration of researchers at Auburn University & the University of Strathc l^de 
I??? 1 /!) . and it is possible to include the states explicitly in the statis- 



tical equilibrium equations (at the expense of potentially very large 
model atoms). 

The sources of atomic data are varied, and in many cases mul- 
tiple data sets for a given ion are available. In some cases these mul- 
tiple data sets represent an evolution in data quality and/or quantity, 
while in other cases they represent different sources and/or compu- 
tational methods. Comparisons of models calculated with different 
data sets and atomic models potentially provide insights into the 
sensitivity of model results to the adopted atomic atoms and mod- 
els (although such calculations have yet to be undertaken for SN). 
Oscillator strengths fo r CNO eleme nts were originally taken from 
iNussbaumer & Storevl dl983l 1 1984) . These authors also provide 
transition probabilities to states in the ion continuum . The largest 
source of oscillator data is from lKuruczl d2009l l2010); its principal 
advantage over many other sources (e.g., Opacity Project) is that LS 
coupling is not assumed. More recently, non-LS osci llator strengths 
have become available through the Iron Project Mummer et al.l 
1993), and wo rk don e by the atomic-data group at Ohio State Uni- 
versity jNaharJl20ld). Other important sourc es of radiative da ta for 
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Fe include Bec ker & Butlerl Jl992L Il995alibh . iNaharl Jl995t) . En- 
ergy levels have generally been obtained from NIST. Collisional 
data is sparse, particularly for states far from the ground state. 
The principal source for collisional data a mong low lying states 
for a variety of s pecies is the tabul a tion bylMendozal Il983l); othe r 
sources include iBerrington et al. dl985J), [ Lennon et alj Jl98 5|). 
iLennon & Burkel Jl994l) Ishine & LinskjT l l 19741) . iTavall Jl997al lbl). 
Zhang & Pradhan ( fl995al ; Il995bl: Il997l). Ph otoionization data is 
taken from the Opaci ty Project " dSeatonlll987l) and the Iron Project 
faummeret"aT] [l993). Unfortunately Ni and Co photoionization 
data is generally unavailable, and we have utilized crude approx- 
im ations. Charge exchange cross-sections are from the tabulation 
bv lKingdon & Ferlanj Il99fj|) . 



To step in time, we adopt a fixed logarithmic increment, 
i.e., At/t = 0.1. Hence, to cover from the initial time of the 
Woosley model of 0.27 d, to the final time of 20.8 d, we computed 
log(20.8/0.27)/ log(l.l) ~46 consecutive timesteps. Each con- 
secutive computation taking on average 1-2 days, the whole time 
sequence takes about 2-3 monthsQ It thus represents a consider- 

2 It is tempting to try to step in time with larger increments (e.g., 20 to 
30% of the current time) but tests we have performed show that the results 
then differ slightly from those obtained with a choice of 10%, or less, and 
are thus not converged. Furthermore, the convergence speed of the code 
is also much poorer since the changes in radiation and gas properties are 
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able computational effort. For the first timestep in the sequence, 
we solve the radiative-transfer problem by fixing the temperature 
everywhere to what it was in the hydrodynamical input. This is 
suitable at depth, but the two solutions differ sizably in the surface 
layers where the assumption of LTE made by KEPLER is invalid. 
Hence, for the second timestep in the sequence, when we let the 
temperature vary and relax as needed, we observe a sizable jump 
in surface properties. Below, when describing our results, we thus 
exclude that original relaxation step, and start our discussion with 
the timestep at 0.3 d, the second epoch in our sequence. 

The spatial grid used for the radiative-transfer solution is set to 
have a constant increment in optical-depth, with typically 5 points 
per optical-depth decade, and a maximum of 140 points. This is 
much less than the original 523 data points of the hydrodynamical- 
input grid. However, our grid is much finer for the regions where 
radiation decouples from matter, the transition from Eddington fac- 
tors of 1/3 to 1 occurs over a few tens of points, while in the hydro 
input this transition takes place over merely a few points. Further- 
more, as the ejecta recombine, our grid adjusts to track the recom- 
bination fronts (across which the optical depth varies steeply). In 
SN1987A, we have three main fronts, associated with the recom- 
bination to He II, He I, and H I. Within these fronts, the stiff and 
complicated coupling between the populations and the temperature 
is one of the main factors slowing down code convergence^ How- 
ever, regions of sharp discontinuities in the hydrodynamical input 
(in composition in particular) are not resolved well and our grid 
then tends to smooth the corresponding variations. Given that these 
gradients are likely overestimated by the ID nature of the hydrody- 
namical model, this may not be such an issue. 



4 EJECTA EVOLUTION 

The simulations we present here were performed over a relatively 
short time, starting at 0.27 d after breakout and ending at 20.8 d. 
While it may seem a short time compared to the bolometric evo- 
lution (with, for example, a light-curve peak at ~80d), this cor- 
responds to an increase in radius by a factor of 77 for all mass 
shells. Hence, over the 46 timesteps computed, the SN ejecta ex- 
panded enormously, and the associated changes at the photosphere, 
as well as above and underneath it, are large. Throughout the dis- 
cussion, we define the SN photosphere as the ejecta location where 
the inward-integrated electron-scattering optical depth is 2/3. In the 
following section, we first present the evolution of the ejecta prop- 
erties, before presenting synthetic spectra and light curves. 

To illustrate the drastic evolution of the ejecta conditions over 
this ~21-day span, we show in Fig. |2]the radius, the velocity, the 
temperature, and the electron density at the photosphere (i? p hot, 
Vphot, Tphot, and n e , P hot) and as a function of time after shock 
breakout. These quantities have been scaled to fit within the same 
plotting range. They are also given in Table [3] together with the 
base ejecta optical depth (taking electron-scattering as the only 
opacity source for the purpose of that table) and the emergent bolo- 
metric luminosity. 

much larger. For computing the long-term evolution of a SN ejecta, one can 
expedite the simulation by starting the time sequence at a few days after 
explosion. 

3 We encounter this severe problem for the present model for SN1987A 
as well as in Type lb SN ejecta models, while for Type II-Plateau SN 
ejecta conditions, recombination fronts develop but never become very 
steep (Dessart & Hillier 2010, in prep). 



Table 3. Summary of properties for the non-LTE time-dependent calcu- 
lations presented in this work. For each epoch, we give the total ejecta 
electron-scattering optical depth (in 10 3 ), the photospheric radius (in 
10 15 cm), the photospheric velocity (in 10 4 km s _1 ), the photospheric 
temperature (in 10 4 K), as well as the emergent bolometric luminosity (in 
10 7 L Q ). 
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During this 21-day timespan the photospheric radius in- 
creases monotonically by over a factor of 10 from 6.0xl0 13 to 
8.1xl0 14 cm, while the photospheric velocity drops from 23 000 
to 4500 km s" 1 . The decrease in photospheric velocity testifies for 
the recession of the photosphere in mass. The photospheric tem- 
perature drops precipitously from 38 100 K due to the fast expan- 
sion, but plateaus after 5-10 days at 5200 K, which corresponds 
to the hydrogen-recombination temperature - it will then remain at 
this value until the nebular phase. Note that our adoption of a floor 
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Figure 3. Left: Temperature distribution as a function of ejecta velocity and time (a color-coding is used to differentiate the epochs). We overplot with a 
dot the corresponding location of the photosphere. From left to right, the two steep drops in the temperature curve at the last time in the sequence (red 
curve) correspond to the He I and H I recombination fronts. Right: Same as left, but now for the quantity R X T, which should be a constant of time for 
an adiabatically-expanding ejecta (we use a logarithmic abscissa scale for better visibility). Deviations from the initial curve stem from energy gain due to 
radioactive decay at depth (below 2000 km s _1 ) and from energy losses (radiation leakage) in the photosphere region and above. The straight line we obtain 
at large velocities and late times results from the setting of a floor temperature of 5000 K up to 0.45 d, then of 3000 K up to 2.8 d, and finally of 2500 K for all 
later times. 




5000 10000 15000 20000 1000 10000 

Velocity [km s" 1 ] Velocity [km s" 1 ] 



Figure 4. Left: Same as for Fig. [3] but now for the electron density. Notice the steady decrease of the electron density with time t, following a general 1 /t 3 
evolution, but steepening at late times and at intermediate velocities as hydrogen recombines in those regions. Right: Same as left, but now for the quantity 
R 3 X n c , a quantity that is a constant of time under fixed ionization. While the curves overlap at early times, several dips develop at later times, associated 
with the recombination of hydrogen (one recombination front) as well as of helium (two recombination fronts). 



temperature (originally 5000 K, and then decreased to 3000 K after 
0.45 d, and down to 2500 K after 2.8 d), used to avoid exponential 
overflow for evaluations involving very high ionization species, is 
well below the temperature of the photosphere at any time, and does 
not affect its properties. 

The electron density at the photosphere varies by about a 
factor of ten during the 21-day span (green curve in Fig. O. 
This is because the optical depth at the photosphere is fixed by 
f_ nndr = 2/3, where k is the continuum mass-absorption 

JK phot 

coefficient and p is the mass density. Its position adjusts as the 
ejecta expand, scaling as the inverse of the photospheric radius. 
In the meantime, the total ejecta optical depth has decreased 
from 2.3 xlO 7 to 3700 (we quote the values for the electron- 
scattering optical depth, which differs little at such times and in 



such hydrogen-rich ejecta from that associated with the Rosseland- 
mean opacity). 

The photosphere location lies, in velocity space, above 
4500 km s _1 at all times, and thus the composition at the photo- 
sphere has not changed through the 21-day evolution (see Fig. [TJ. 
At the photosphere, we have the following mass fractions : Xn = 
0.478, Xhc = 0.517, X c = 4.2 x 10" 5 , X N = 2.5 x 10" 3 , 
Xo = 4.9 x 10~ 4 , X Fc = 7.9 x 10~ 4 , X m = 3 x 10~ 5 (note that 
the Nickel is primarily in its stable form, the corresponding Ni 
abundance being initially 3 x 10~ 8 ). 

Despite the relatively small time span of 21 days of our time 
sequence, it covers a fast and very rich evolution for SN1987A. 
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Figure 5. Same as Fig. [5] but now for the level of ionization of hydrogen (top left), helium (top right), oxygen (bottom left), and iron (bottom right). An 
ionization level n, where n is a positive real number, for a species X means that this species X is found at that location primarily in its n th -time ionized form 
(i.e., the ionization level is equal to ^ iX l +/ Jj. X l +) (see text for discussion) 



4.1 Ejecta temperature 

The evolution of the ejecta temperature distribution echoes what 
was shown above for the photospheric temperature. In the inner- 
most part of the ejecta, below ~ 2000 km s , the temperature drops 
from about 10 6 to 5x 10 4 K, while above that, it is on the order of 
1.5 x 10 4 K or less. This is a huge change from the early very hot 
ejecta conditions we started with (Fig. [3j dots give the location of 
the photosphere at each time). Notice the steepening of the tem- 
perature distribution just below the photosphere, associated with 
the hydrogen recombination front. Two other fronts also develop, 
a strong one for the He I-He II transition and a weaker one for the 
He II-He III transition, both of which occur at large optical depths. 

To reveal the ingredients that control the temperature we inte- 
grate Eq.[3]over frequency to give 

1 D(r 4 J) 1 d(r 2 H) [°° . T , J 

Since at depth, J m B, J 2> H, and the integral on the right- 
hand side is approximately zero (ignoring radioactive heating) we 
have that r 4 J is constant, and hence r x T is a constant of time. 
The same results follow from Eq. [2] if we treat the radiation as a 
photon-gas with e = aT 4 / p and P = aT 4 /3. 

Deviations from that constant r x T curve will reflect depar- 
tures from adiabaticity, associated here with the energy decay (a 



gain in regions where 56 Ni is present) and radiation losses (near 
the photosphere where escaping photons extract energy from the 
thermal pool). But regions lacking 56 Ni and situated at large opti- 
cal depth should evolve at constant r x T. 

In the right panel of Fig. [3] we show the time evolution of 
that quantity r x T obtained in our simulations, plotted for all 
ejecta mass shells as a function of time (we use a logarithmic scale 
for the abscissa for visibility). At early times, all the curves over- 
lap, but within a few days, a pronounced "bump" forms below 
2000 km s _1 , an excess which is caused by the energy contribu- 
tion from radioactive decay. Importantly, we note that it does not 
extend beyond 2000 km s _1 , implying that the locally-deposited 
energy from decay has not had time to diffuse out, even at 20.8 
days after explosion. As the photosphere is located at (or above) 
4500 km s - , radioactive decay has no effect on our synthetic spec- 
tra during these initial 20.8. 

At large velocities, the radiation that escapes from the photo- 
sphere makes the curves dip. Note that the deficit extends below 
the photosphere, in regions that "feel" the energy leakage taking 
place immediately above. While Fig. [3] highlights some key prop- 
erties governing the energy budget of SN ejecta in general, and 
of SN1987A in particular, it also serves as an important check for 
the numerics. Indeed, in regions where the 56 Ni mass fraction is 
small or negligible, and located well below the photosphere, the 
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expansion is adiabatic and we find that the temperature scales as 
the inverse of the radius, as it should. 

4.2 Ejecta ionization 

We now turn to the evolution of the electron-density distribution 
as a function of time (Fig. |4). Due to the homologous expansion, 
the mass-continuity equation leads to a 1/t 3 evolution of the mass 
density at a fixed velocity. Furthermore, if the evolution is at con- 
stant ionization, the quantity r 3 x n a is a constant of time. We plot 
this quantity r 3 x n c in the right panel of Fig. [4] While the curves 
overlap at early times, there is a systematic and growing deviation 
at large velocity, starting just below the photospheric velocity at 
each time. This deviation starts when the conditions of full ion- 
ization cease to prevail. It is caused by hydrogen recombination, 
which severely depletes the density of free electrons. At depth, a 
similar but more modest deviation occurs, and is associated with 
the recombination of helium. 

While hydrogen and to a lesser extent helium recombination 
impact the electron density visibly, all species present in the ejecta, 
no matter how relatively under-abundant, will recombine from their 
high ionization states at the time of breakout to a neutral or once- 
ionized state after a few days/weeks. We illustrate this systematic 
recombination in Fig.|5]for hydrogen (top left), helium (top right), 
oxygen (bottom left), and iron (bottom right). As before, a color 
coding is used to differentiate the epochs - recall that the time incre- 
ment is constant in the log (it increases by ~10% at each time step). 
The key characteristics that apply to all is that the ionization is high 
at low velocity, reaches a minimum around the photosphere, but in- 
creases again at high velocity (in particular for the dominant species 
hydrogen and helium). This behavior reflects the ionization freeze- 
out in fast-expanding low-density SN ejecta, and stems from in- 
corporating time- dependent terms in the st atistical-equilibrium and 
energy equations jDessart & Hiilierll20080 . This ionization freeze- 
out cannot prevent recombination altogether, but the effect is large 
enough to dramatically alter the strength of numerous strong lines, 
such as H I Balmer lines or Na ID (see jj6}. 

Initially, the material at depth is highly-ionized for all species, 
i.e., H + , He 2+ , 5+ , Fe 7+ , while at the photosphere we have 
H + , He 2+ , 4+ , Fe 5 ' 5+ . These are highly-ionized conditions that 
were not seen at the photosphere of SN1987A because observations 
started "too late", at one day after explosion, when the ions that pre- 
vail at the photosphere are: H + , He + , 2+ , Fe 3+ . At the end of the 
simulation at 20.8 d, the material at depth is still hot and optically- 
thick, and we find H+, He + -He 2+ , + ~0 4+ , Fe + -Fe 5+ , while 
at the photosphere the ions are now close to their neutral state: 
H 0+ , He 0+ , O 0+ , Fe + . This strong contrast in ionization between 
the optically-thick hot inner layers and the optically-thin cool outer 
layers is what forces us to retain a broad range of ionization stages 
for all important ions. In our current brute-force approach, any ion 
important at some depth must be included at all depths. 

The photosphere location tends to track the ionization mini- 
mum, bracketed between the fully ionized hotter layers at large op- 
tical depth and the cool outer optically-thin regions subject to ion- 
ization freeze-out. So, in the new approach followed, we retrieve 
the properties identified in that former work, but now within the 
context of a calculation based on a physical hydrodynamical input 
of the explosion and pre-SN progenitor evolution, and without any 
prescribed boundary conditions for the flux. The evolution we see 
here further implies that events like SN1987A experience a faster 
change in ejecta ionization within their first day, motivating the col- 
lection of multiple spectra within the first day/night after explosion. 



5 SYNTHETIC SPECTRA AND BOLOMETRIC 
LUMINOSITY 

In this section, we present the radiative properties of SN1987A, 
focusing exclusively on the results from our non-LTE time- 
dependent radiative-transfer simulations. We s cale our synthetic 
spectra adopting a distance of 50 kpc to the LMC dRibas et al-120021 : 
lDall'Oraetalj2004 IClement et ai]|2008l ; IPietrzvfiski etall2 009 ) . 
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Figure 6. Montage of synthetic spectra between 0.30 d and 4d. For better visibility, we interpolate in time between those epochs that were computed, and show the spectra at the following times: 0.30, 0.40, 0.53, 
0.71, 0.95, 1.26, 1.69, 2.25, 3.00, and 4.00 d. We used a constant increment in the log to better resolve the fast evolution at early times. No reddening is applied, and the only scaling is for the distance, taken at 
50kpc. 




Figure 7. Same as Fig. [6] but now for later times. The split in time between the previous figure and this one corresponds to the epoch when the He I 5879A line vanishes. This transition occurs prior to the appearance 
of NalD. We use a constant time increment, and we plot synthetic spectra at the following times: 4.00, 5.87, 7.73, 9.60, 1 1.47, 13.33, 15.20, 17.07, 18.93, and 20.8 d. No reddening is applied, and the only scaling is 
for the distance, taken at 50 kpc. Notice how Ha strengthens with time despite the declining Lyman/B aimer continua and the cooling of the photospheric layers down to the hydrogen-recombination temperature. 
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5.1 Spectroscopic evolution and sources of line blanketing 

As our time sequence starts at 0.27 d after explosion, we compute 
non-LTE time-dependent spectra and light curves starting at epochs 
that were not observed (i.e., before day 1) and extending to 20.8 d. 
We show our synthetic spectra for the early/fast evolution in Fig. [6] 
(day 0.3 to 4; we omit the first date in the sequence, which was 
computed at fixed temperature, as given by the hydrodynamical in- 
put) and the subsequent/slow evolution in Fig. [7] (day 4 to 20.8). 
Given the BSG nature of the progenitor, and in particular the mod- 
est progenitor-star radius of ~50Rq at explosion, the initial fast 
expansion of the shocked outer layers of the progenitor star in- 
duces a dramatic cooling. This is the cause of the much faster spec- 
tral evolution of SN1987A compared to Type II-P SNe which have 
RSG progenitor stars with surface radii of 500-1500 Rq. 

Our simulations give a photosphere whose temperature de- 
creases from 38 100 K at 0.3 d to ~ 1 1 000 K at 2.0 d. This is then 
followed by a slow decrease to >5000K over a number of days. 
As expected, and as shown in Fig. [6] we obtain synthetic spectra 
whose peak-flux shifts very rapidly during these initial two days, 
from around 300A to 3000A. Over that time, the UV flux drops by 
about two orders of magnitude while the optical flux rises by a fac- 
tor of a few. The UV flux is initially highly wavelength-dependent, 
caused by strong line blanketing from multiple ions. In contrast, 
the optical suffers little from line blanketing at early times and the 
corresponding flux is essentially featureless. 

Using the photosphere as a reference location, we show in 
Fig. [8] the spectral regions where important line blocking occurs 
at three epochs, differentiating the contributing ions with a color 
coding. At 0.33 d, line blocking appears at short wavelengths (200- 
2000A) and is associated with the highly-ionized ions that are 
present at the photosphere at those early times, i.e., O III-V, Si IV, 
S IV, Fe V- VII, Co VI- VII, Ni VI- VII (top panel of Fig. [8}. Hence, 
considerable blanketing takes place at such times and under highly- 
ionized conditions, although this takes place in the UV and is rarely 
seen. Lines in the optical region are extremely weak (in part be- 
cause of the very large photospheric velocity at those early times), 
and are associated with H I (Balmer lines), He I and He II. 

At 3.05 d (middle panel of Fig. [8}, the spectral regions with 
strong line blocking are still found shortward of 3000A, but the 
associated ions are of a lower ionization state (e.g., Oil, Si II- III, 
S II- III, Ti III, Fe II- IV, Co II- III). Ni II and Ni III would also con- 
tribute but have not been included in this calculation. In the optical, 
a small number of lines (relative to the forest of metal lines) asso- 
ciated with H I and He I make their conspicuous appearance. He II 
4687A is gone, and lines of once-ionized CNO elements are now 
present. 

Beyond 3-4 days, the SED continues its dramatic fading in 
the UV region (see bottom panel of Fig.[8]for an example at 9.8 d), 
which decreases by a total of about five orders of magnitude from 
0.3 d to 20 d. In contrast, past 3^t days and after the disappearance 
of He I 5 877 A, the optical range evolves slowly with a strength- 
ening of the hydrogen Balmer lines, Call, and metals of low- 
ionization state (Sill, Till, Fell, Coll). Our simulations do not 
include Ni II-IV in the model atom, and have zero abundance for 
Scandium and Barium. Hence, no microphysical processes associ- 
ated with these species/ions are predicted by our computations. 

We wish to emphasize that at all times line blanketing severely 
depletes any overlapping flux. It achieves this very effectively in the 
UV, in association with ions that may have a high or low ionization 
state. Line blanketing is thus not born when conditions are cool; it 
is instead the case that as conditions cool, the spectral region where 
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Figure 8. Top: Illustration of the wavelength distribution of all lines having 
a Sobolev optical-depth Ti; nc greater than 0.01 at the SN photosphere and 
0.33 days after the explosion. We compute Ti; nc in the co-moving frame 
using the Sobolev approximation. The abscissa range is 100-2000A. This 
illustration does not intend to show the effective opacity felt by photons, 
but rather how line opacity is distributed with wavelength and between ions 
(we overplot all line contributions rather than add them). Middle: Same 
as top, but now for the SN at 3.05 d. The abscissa range is 100-7000A. 
Bottom: Same as top, but now for the SN at 9.6 d. The abscissa range is 
2000-10000A. 
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Figure 9. Illustration of the evolution of the photospheric velocity (black) 
in comparison with the velocity measured at maximum absorption in the 
P-Cygni profile troughs we associate with Ho, H/3, Hy, and Fell 5169A. 
Measurements are done on synthetic spectra. At times prior to 3-4 days, 
line overlap tends to affect very weak lines. Hence, the measurement may 
be a hybrid of multiple contributions at earlier times. However, notice how 
well Fe II 5169A, once it has become strong enough beyond a few days after 
explosion, traces closely the photospheric velocity. Notice also how all lines 
u nderestimate the p hotospheric velocity at early times, a property discussed 
in lDessart & HillieJ J20053) . 

line blanketing occurs shifts towards (observed) optical spectral re- 
gions. One must also recall that as this takes place, the peak of 
the flux distribution shifts to longer wavelengths, so that there is 
less and less flux to block. For example, at 20 d, the line blanketing 
caused by Fe II is enormous at UV wavelengths, but there is hardly 
any photons injected at the thermalization depth in the UV at such 
times. 

5.2 Ejecta kinematics from line profile morphology 

Inferring the photospheric velocity of the SN ejecta from spectro- 
scopic observations at a given epoch is important as it enters the 
determination of the total kinetic energy of the ejecta. It is also 
a key component for the determination of dis tances to Type II 
SNe using the expanding-photosphere method dKirshner & Kwanl 
ll974l ; lEastman et al.|[T996l : iDessart & HilliertkoOSd). the s pectral- 
fitting expanding-atmosp here method ( Mitchell et al. 2002), or the 
standard-candle method drlamuv & Pintoll2002l) . In Fig. [9l and us- 
ing the synthetic spectra for our present time sequence, we show the 
Doppler velocity corresponding to the location of maximum profile 
absorption for several optical lines (Ha: Magenta; H/3: Blue; H7: 
Green; Fe II 5 169A: Red). We find that at early times all lines tend 
to show a maximum absorption at a Doppler velocity that under- 
estimates the photosp heric velocity (shown in b lack), for a reason 
discussed in detail in Dessart & Hillier] d2005al) . However, once a 
given line strengthens, that Doppler velocity then overestimates the 
photospheric velocity, by an amount that varies from being large 
for Ha, to modest for H/3 and H7. For Fell 5 169 A, this Doppler 
velocity underestimates the theoretical photospheric velocity at all 
times prior to ~15 d, and modestly overestimates it afterwards. 

In a separate test calculation, not shown here, we evolved 
the same hydrodynamical model but with a three times larger iron 
abundance. In that case, the Fell 5169A-line was stronger at all 
times and the corresponding velocity location of maximum profile 
absorption matched more closely the photospheric velocity. The ac- 
curacy of using the Fell 5169A line for inferring the photospheric 
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Figure 10. Bolometric light curve for SN1987A based on our non-LTE 
time-dependent radiative-transfer calculations with CMFGEN (black; dots 
refer to the epochs computed). We also overplot the prediction, for the same 
ejecta, obtained with V1D I Dessart et al. 2009), which agrees to within 10- 
20%. Both curves ag ree qualitatively and quantitatively with the results of 
iBlinnikov etaD feOOfl see their Fig. 8). 

velocity is thus conditioned by the environmental metallicity, and 
may therefore vary with redshift or galaxy hosts. 

5.3 Synthetic bolometric light curve 

To conclude this section on synthetic radiative properties, we 
present the bolometric light curve for our time sequence in Fig.llQI 
Starting at 0.3 d, we miss the breakout phase, but capture the steep 
fall-off that immediately follows it until the luminosity levels off 
after about 7 days. Beyond that, we then see a slow increase in 
bolometric luminosity. 

In Fig. [TO] we overplot the prediction of the one -dimensional 
radiati on-hydrodynamics V1D code, as described in lDessart et al.l 
(2009). This computation is done assuming gray transport (i.e., 
one group), a mean opacity that neglects line contributions, and 
flux-limited diffusion. The bolometric luminosity obtained with 
this much coarser approach agrees to within ~20%. Furthermore, 
our non-LTE time-dependent line-blanketed treatment, shows an 
evolution that agrees with the inferred observed bolometric light 
curve as well a s the r adiation-hydrodynamics results obtained by 
IBlinnikov et"aH d200(j| . see their Fig. 8 for comparison with both 
their models and the inferred bolometric luminosity). 

The luminosity plateau and subsequent increase after ~7 d 
occurs at a time when the 56 Ni energy which, in this model, we 
deposit at depths below ~ 2000 km s _1 (there is a buffer of about 
10 Mq between the photosphere and the 56 Ni-rich regions at that 
time), has not had time to diffuse out. The "instantaneous" diffusion 
time for these inner and very dense layers is indeed orders of mag- 
nitude larger than the age of the SN. Hence, this plateau and subse- 
quent re-brightening is not caused/triggered by energy deposition 
from unstable isotopes, but results instead from the rate at which 
the photospheric radius and temperature happen to evolve in the 
shock-heated envelope. Seven days after explosion corresponds to 
the time in our simulations when the photospheric temperature be- 
comes essentially constant at >5000 K, which corresponds to con- 
ditions under which hydrogen recombines to its neutral state. At 
the same time, the photospheric radius continues to grow and the 
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Figure 11. Observed (filled circles) and theoretical (lines) U BVRI light curves as a function of time since explosion. We restrict the time baseline to 21 days 
to focus on that period over which we compute our time sequence. The source of observational data is composite, coming from CTIO obs ervations an d IAU 
circulars for earlier times. The synthetic photometry was computed by integrating our synthetic spectra over the Landolt filter bandpasses (Landolt 1992]), and 
applying a reddening of E(B — V) = 0.15mag and scaling for a distance of 50 kpc. We plot crosses along the synthetic [/-band light curve to indicate the 
epochs at which our calculations are carried out. Notice the excellent match between theory and observations for the BVRI bands. For the U -band, our values 
lie above the observed ones at all times (see text for discussion). 



flux has to go up (see Fig.[2]for an illustration of those trajectories). 
The photospheric velocity at that time is ~9000 km s _1 and there is 
no dramatic change in density gradient at that location. This point 
of inflection in the bolometric light curve would likely be shifted 
for progenitor stars with different initial radii or for explosions of 
differing magnitude. 



6 COMPARISON TO OBSERVATIONS 

Having discussed the radiative properties of the synthetic spectra 
of our SN1987A time sequence, we now confront the se to observa- 
tions. We use the observational dataset presented in Philli ps et al.l 
Jl988h . complem ented with early-time photometric data published 
in IAU circulars feunkel et al]|l987l ; iHerald et al.lll987T). and th e 
UV spectra obtained with IUE and presented in lPun et alj Jl995l) . 
We sc ale our synth etic s pectra adopting a dist ance of 50 kpc to the 
LMC llRibas et alJl2002l ; lDall'Ora et alj|2004 IClement et alj2008t 
IPietrzvnski et alj2009l) . When corrected for extinction, as when we 
compare to observed spectra and magnitudes, we emp loy a redden- 
ing E (B — V) = 0.15 and the reddening law of ICardelli et al.l 
( 1989). This seems adequate for obtaining a satisfactory match to 
observed spectra, although this value could be altered by ±0.05 and 
remain suitable. 



6.1 Observed versus synthetic multi-band light curves 

In Fig.QT] we show the observed (filled circles) and synthetic (lines 
and/or crosses; adopting a reddening of E(B — V) = 0.15 and a 
distance of 50 kpc) UBVRI light curves for SN 1987A. The mag- 
nitudes are computed by integrating t he synthetic s pectra over the 
Landolt filter-transmission functions (Landol3 ll992h . We obtain an 
excellent agreement for the optical bands BVRI, but our predic- 
tions for the [/-band magnitude is offset, predicted too bright at 
all times. This discrepancy could be somewhat reduced, although 
not resolved, by choosing a slightly higher reddening of 0.2 mag. 
The (7-band magnitude decrease is also too slow in our simula- 
tions, and levels at a minimum that is ~1-1.5 mag too bright. This 
mismatch may be physical, perhaps reflecting a slight offset in the 
density distribution in the progenitor surface layers. The progenitor 
star model should perhaps have a slightly smaller radius, resulting 
in a faster evolution of the SED to longer, optical, wavelengths. It 
may also be associated with opacity issues, stemming from under- 
abundant but nonetheless important species we presently ignore. 
For example, we are not including Al, which has strong resonance 
lines in the UV, or, e.g., Cr, Ar, Ne which do contribute opacity 
there. We may be missing some low ionization stages, e.g., Nill 
or Fe I, although the U -band magnitude mismatch is present at all 
times, not just at late times when the photospheric conditions are 
cool. Finally, we may be using a too small model atom for some 
species. For example, we include only 115 levels of Fell, which 
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we find yields converged results at early times, but is not quite suf- 
ficient when Fe II is the dominant Fe ion (we find that enhancing the 
number of Fell levels from 115 to 827 at the recombination epoch 
(past day 5) leads to a reduction of the already extremely small flux 
in the 1000-2000 A range, but leaves the spectrum at longer wave- 
lengths unaffected). Recall here that iron has a mass fraction of 
less than 0.001 in our model, while hydrogen is the most abundant 
element and is the main donor of the free-electrons that provide 
the bulk of the opacity. In any case, the mismatch cannot be asso- 
ciated with something sizably inadequate in the progenitor model 
or in the code since we obtain excellent agreement for the optical 
bands. It also intervenes in regions where the flux is orders of mag- 
nitude smaller than in the optical, or orders of magnitude smaller 
than what it used to be at 0.3 d. But this goes to show that even in 
the case of SN1987A, for which we have tight constraints on the 
progenitor star and the explosion properties, obtaining a complete 
agreement is a challenge. 

6.2 Observed versus synthetic spectra 

We present in the left panel of Fig. [12] a montage of all observed 
UV and optical spectra of SN1987A up to 20 days after explo- 
sion. These are stacked up vertically and shown on a logarithmic 
scale for visibility. A color coding differentiates the epochs. In the 
right panel, we show a montage of synthetic spectra (reddened us- 
ing E(B — V) = 0.15 and scaled for a distance of 50kpc), inter- 
polated in time from our time sequence so that the displayed dates 
match each observational epoch in the left panel. Importantly, the 
same vertical shift is applied to both observed and synthetic spec- 
tra for the same epoch. Overall, the synthetic spectra shown in the 
right panel mirror with high fidelity the observed spectra shown in 
the left panel. 

For the first few epochs, echoing what was seen already in 
Fig. QT| the synthetic fluxes are somewhat larger than observed. 
The color temperature is slightly too large, as is the ionization state 
of the photosphere, so that all lines appear weaker above the con- 
tinuum. If physical, this mismatch could reflect the fact that the 
progenitor star used as an input in this simulation has a slightly too 
large initial radius. A more compact progenitor star would induce 
a faster cooling and a faster UV fading. 

After a few days, the agreement between observations and the- 
ory is excellent, in particular in the optical where both continuum 
and lines fit to within ~10% at all times. There are a few excep- 
tions to this, but Ha is one line that embodies this agreement, as it 
is reproduced in both strength and width. In the UV, the synthetic 
flux overestimates the UV flux by a few tens of percent at most 
times. This echoes the flux mismatch in the U band although the 
origin may not be the same (opacity sources in the UV and around 
3000-4000A differ; see Fig.©. Recall that the UV flux drops by 5 
orders of magnitude between 0.3 and 20.8 d, and that the UV flux 
after a few days is a very small fraction of the total flux. In any 
case, we will investigate this issue in the future, but here, given the 
complete lack of free parameters after we start the time sequence, 
this agreement is impressive. 
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Figure 12. Left; Montage of observed spectra for SN1987A, including the optical at all times and the UV when available ^Phillips et al . 1988; Pu n et alj|l995t> . For visibility purposes, we apply a scaling between 
epochs, ordering the spectra chronologically starting at the top. A color coding is used to differentiate each date. Right: Same as left, but now showing our synthetic spectra, reddened with E(B — V) = 0.15 and 
scaled for a distance of 50kpc. Since our time-sequence was not computed for the specific observed times, we interpolate in between the computed epochs to make the correspondence exact in time between left and 
right panels. Moreover, we apply the same final scaling, and the same color coding, as in the left panel to facilitate the visibility. Importantly, apart from a small offset for the first epochs, the observed and theoretical 
fluxes agree to within about 10 % in the optical (the agreement is not as good in the UV range), which is quite remarkable given that we started from a hydrodynamical model of the explosion and just let the ejecta 
evolve "naturally" from day 0.27 onwards. 
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We now present detailed fits to individual observed spectra. 
For this exercise, we allow ourselves to choose a synthetic spec- 
trum whose date may be offset by ±0.5 d from that of the observed 
date. This matters at early epochs, when the evolution is fast, and 
for later dates when our photosphere seems to have somewhat too 
high ionization/temperature (the [/-band flux is visibly too high). 
Recall that we cannot tune the model parameters to fit the observa- 
tions since all the ejecta properties were inherited at the start of the 
sequence from the adopted hydrodynamic explosion model. 

The first two spectra obtained for SN1987A on the 24th and 
the 25th of February 1987 show a blue continuum, which repre- 
sents a fossil of the conditions at breakout (top row of Fig.|13ll. We 
see broad and strong B aimer lines as well as the conspicuous pres- 
ence of He I 5 877 A, reproduced here with the helium mass frac- 
tion prevailing at the surface of the B SG progenitor star. Hence, i n 
contrast with the pre v ious results of Eastman & Kirshnerl (1989); 
ISchmutz et all Jl990h ; iHQflicbl jl987t Il988h . it is indeed possible 
to reproduce this line with an abundance typical of a BSG progen- 
itor in the LMC. While it is difficult to identify the specific ingre- 
dients) resolving this former discrepancy, this highlights the need 
for realistic radiative-transfer models for an accurate reproduction 
of SN spectra, accounting for non-LTE effects, line blanketing, etc. 

The fit to observations is fair given the very fast evolution of 
the photospheric conditions at these early times and the strong de- 
pendency on progenitor properties (in particular radius). More fre- 
quent observations over these initial two days would have revealed 
further detail about this fast evolution and provided additional con- 
straints for the modeling. 

By the third observation on the 26th (middle-left panel in 
Fig. [13}, Hel5877A is gone and the first signs of metal lines in 
the optical appear (e.g., FeIl5169A). In our calculations, we still 
predict He 1 5877A (although it is no longer present in our time se- 
quence one day later). Another sign that the ionization is too high 
at such times is the overestimated UV flux. The Ha line profile is 
well fitted, although the emission part is somewhat underestimated. 
Recall that the early evolution is fast and that a very optically-thick 
line like Ha is sensitive to even slight modulations in ionization 
(ionization is changing rapidly at this epoch as hydrogen recom- 
bines in the photospheric layers). Interestingly, the H/3 line is not 
as well fitted as Ha, although its formation process is analogous, 
and this discrepancy lasts for about a week. 

Continuing onwards, the ejecta recombine at the photosphere 
and the UV flux essentially vanishes. We fit the Ha line pro- 
file and its evolution very well, both quantitatively in strength 
and width, but also qualitatively in its shape (i.e., how the flux 
varies from the blue edge to the red edge of the profile). This 
suggests that, together, our calculations and the initial model al- 
low an accurate description of the Ha line-formation process; the 
density distribution and the kinetic energy of the original model 
are well suited. This contrasts with the former discrep ancy in 
Balmer-line strengths encountered by |Schmutzetal.Ul990h . which 
at the time were associated with the potential role of clumping 
(note that no clumping is used here). The good reproduction of 
Balmer-line strengths or NalD supports the adequacy of our com- 
puted non-LTE time-dependent ionization structure. As pointed 
out by Utrobin & Chugai J20051) . and independently confirmed by 
iDessart & Hillied d2008h . the competing recombination and expan- 
sion timescales for the ejecta lead to an ionization freeze-out, par- 
ticularly visible in the optically-thin fast-expanding low-density re- 
gions. We capture this effect by including time-dependent terms 
in the statistical-equilibrium and energy equations; this important 
aspect requires a full non-LTE treatment. 



These findings have been challenged by |Pe et al.l j2009h but 
their tests on SN1987A out to 20 days were performed using ejecta 
conditions that produce a strong Balmer continuum at late times 
(their figur es 8, 9) which is not seen in IUE observations (pre- 
sented by IPun et al.ll 19951 modeled in this paper, and shown, e.g., 
in Fig. 1 124415b . This suggests that in their simulations the photo- 
spheric layers at such times are well above the recombination tem- 
perature for hydrogen and thus copiously ionized. In their model, 
ionization is likely controlled by Lyman/Balmer continuum pho- 
tons, and since the photospheric layers are not recombining, the 
ionization freeze-out we predict cannot take place in their inves- 
tigation. In contrast, time-dependent effects become important for 
determining the ionization structure of the ejecta when there are no 
Lyman/Balmer continuum photons to contribute, which inevitably 
occurs when hydrogen recombines, at around 7000 K. 

After the 1st of March 1987 (Fig. [T4lTT5V the optical flux be- 
comes better fitted while we overestimate the very-faint UV flux. 
The Na I D doublet appears in our models on March 3rd, about a 
day later than in observations. It strengthens in time but matches 
the observations well after a few more days. The fit to the Ha line 
profile is quite spectacular. Recall that we do not treat Scandium 
and Barium and thus do not include their optical features. 



7 DISCUSSION AND CONCLUSIONS 

We have presented a new approach for radiative-transfer modeling 
of SN ejecta, retaining the key assets of our former approaches. 
By treating non-LTE explicitly and incorporating time-dependent 
terms in the statistical-equilibrium, energy, and moments of the ra- 
diative transfer equations, we improve the physical consistency of 
our computations. In particular, the simulations are performed on 
the entire ejecta, starting from a given hydrodynamical input at a 
given post-explosion time, making this an initial-value problem, 
while spatial boundaries are zero-flux at the inner edge and free- 
streaming at the outer edge. As diffusion at the inner boundary does 
not need to hold we can evolve the ejecta from its photospheric 
to its nebular phase. Importantly, integrating our synthetic spec- 
tra over filter-transmission functions, we compute non-LTE light- 
curves that account for the explicit role of line-blanketing, without 
recourse to the standard expansion opacities and the assumption 
of LT E for the gas, as is typically used (see, e.g.. lBlinnikov et aU 
l200d) . 

Delaying a presentation of all the technical details of this new 
approach (Hillier & Dessart 2010, in preparation), we illustrate in 
this paper these modeling improvements with results for SN1987A, 
using as our starting conditions the ejecta composition and struc- 
ture of the hydrodynamical input model "lml8a7Ad" of Woosley 
(priv. comm.). We evolved the ejecta gas and radiation from 0.27 
to 20.8 d after explosion, corresponding to a spatial expansion of 
a factor of 77, and considerable cooling. At the photosphere, the 
ejecta conditions evolve from fully-ionized to once-ionized and 
neutral. During the 21 day time sequence, the photosphere recedes 
in mass, and thus in velocity, while its radius increases steadily. 
After a few days, the photospheric location is primarily set by the 
location of the hydrogen ionization front. Our computed temper- 
ature and electron-density evolution reproduces predictions for a 
radiation-dominated homologously-expanding gas, modulo the ef- 
fect of cooling at the photosphere and heating from radioactive 
decay in the inner ejecta. At the last time in the simulation, the 
photosphere is at 4500 km s _1 , and is not influenced, either di- 
rectly or indirectly, by energy deposition from radioactive decay 
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Figure 13. Comparison between observed (black) spectra for SN1987A and our synthetic fits (red; reddened with E(B — V) = 0.15 and scaled for a distance 
of 50kpc) from our non-LTE time-dependent sequence, covering the epochs from the 24th of February until the first of March 1987. A small absolute scaling 
of at most 30% is applied to the synthetic spectra to bring them down to the level of the observations (see text for discussion). 



occurring below 2000 km s -1 in our model. The surface conditions 
computed here and prior to 20 d are thus not affected by unstable 
nuclei. Since the SN1987A model ejecta we employ are homoge- 
neous above 4500 km s -1 , the changes we observe in the spectrum 
over this 21 -day time span are conditioned by modulations in ion- 
ization rather than composition. 

Adopting a distance of 50kpc and a reddening E(B — V) = 
0.15, we obtain very good agreement between synthetic and ob- 
served light curves in the optical. Our £/-band synthetic magni- 



tudes are, however, too bright by a fraction of a magnitude at early 
times, and this discrepancy slowly grows as the flux in the blue 
ebbs. Our synthetic spectra computed for the epoch 0.3-1 d cannot 
be compared to observations, which started one day after explo- 
sion. At such early times, we find an SED with a peak flux in the 
far-UV, where strong line blanketing occurs due to 3-4-times ion- 
ized metal species (O, Fe, Ni), while the nearly-featureless optical 
range shows weak lines associated with Hi, Hell and He I. After 
day 1 and until the 14th of March 1987, we compare our predic- 
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Figure 14. Same as Fig. |T3] but showing dates between the 2nd and the 8th of March 1987. The reddened and distance-scaled synthetic spectra are within a 
few percent at most of the observations. 



tions to observations. In general, the flux agreement in the optical 
is at the 10% level (in an absolute sense, i.e., without any normal- 
izati on), while in the UV range it is within 30-50%. In contrast 
with lEastman & KirshneJ fl989), we reproduce the strength and 
morphology of He I lines at early times, under a photospheric com- 
position that is compatible with a BSG progenitor star. Through- 
out the time span of our sequence, we reproduce well the multi- 
band light curve, the overall SED, and in particular the strength 
and shape of Ha. Since radioactive decay does not influence the 
computed photospheric conditions in any way at such times, and 



given the relatively good fit to observations (in particular if we fo- 
cus on the bulk of the emergent luminosity, which falls at opti- 
cal wavelengths), it does not play a pivotal role for understanding 
the radiative properties of SN1987A up to 20 d, in contradiction 
with the proposition of iMitchell et al.l bOOll |2003) . In contrast, 
our good reproduction of the Balmer line profiles at all times, and 
in particular when there are no longer any photons in the Lyman 
and Balmer continua, gives strong support to our proper modeling 
of the time-dependent ionization st ructure of the SN1987A e jecta, 
an important issue emphasized by lUtrobin & Chugail d2005n and 
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Figure 15. Same as Fig. [T3] but showing dates between the 9th and the 14th of March 1987. Some of the prominent discrepancies are caused by missing 
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iDessart & Hillierl d2008h . be et all d2009h suggest that such time- 
dependent effects are much weaker, but their presentation for the 
recombination phase of SN1987A is unconvincing since they use 
models that are fully-ionized and UV bright (in contradiction with 
observations). Furthermore, despite the one-dimensionality of our 
approach, we achieve very competitive fits to observations. One 
may argue the unresolved problems we encounter are caused by 
those neglected multi-dimensional effects, but we surmise this de- 



parture from sphericity would also alter line profile shapes etc., 
which are fairly well reproduced (see, e.g., Ha). 

We predict the fastest spectral evolution in SN1987A occurred 
prior to day 1. From 0.3d to Id, most of the flux was emitted 
in the range 300-3000A which could have been captured by the 
IUE satellite. Prior to 0.27 d, the photosphere would have shone 
at shorter wavelengths, but the UV range would have revealed the 
long-wavelength tail of that SED, testifying for the hotter condi- 
tions at earlier post-breakout times. This suggests that obtaining 
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multiple observations prior to day 1, across as broad a wavelength 
range as possible, would have captured the phase of fastest evolu- 
tion of SN1987A. While the UV shows the largest changes, multi- 
ple spectra, rather than a unique spectrum, of SN1987A during both 
the first and the second night of its evolution would have been valu- 
able0This applies to some extent to Type II-P SN as well, although 
their evolution is slower. In Type I SNe, the SED likely evolves too 
fast after breakout to allow us to capture it, but then, the first few 
days after explosion would be very useful for constraining the sur- 
face composition of the progenitor. Overall, observing frequently at 
early times would provide important constraints that are generally 
missing in current observations. 

The work presented here is an essential benchmarking of our 
code since SN1987A is one of the best understood and best ob- 
served supernova. As it places numerous and tight observational 
constraints, it offers a valuable alternative, although not a re- 
placement, to benchmarking against other codes. The combined 
approach of using non-LTE, time dependence, and line blanketing 
seems very promising. The future from here is to systematize such 
investigations by evolving a wide range of ejecta for all supernova 
types and compare results with multi-epoch multi-wavelength ob- 
servations. The aim is then to use such models and observations 
to make quantitative inferences on the progenitor properties, pre- 
SN star evolution, and the explosion mechanism. In the context of 
SN1987A, we now need to repeat the present exploration with a va- 
riety of progenitor stars and explosion properties, in order to delin- 
eate the systematics associated with core-collapse SN explosions of 
BSG stars. In parallel, it would be valuable to gather additional and 
high quality observational data for similar Type II-peculiar events 
like SN1987A. 
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